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ABSTRACT 

We generalize the classic Bonnor-Ebert stability analysis of pressure-truncated, 
self-gravitating gas spheres, to include clouds with arbitrary equations of state. A 
virial-theorem analysis is also used to incorporate mean magnetic fields into such 
structures. The results are applied to giant molecular clouds (GMCs), and to 
individual dense cores, with an eye to accounting for recent observations of the internal 
velocity-dispersion profiles of the cores in particular. We argue that GMCs and massive 
cores are at or near their critical mass, and that in such a case the size-linewidth 
and mass-radius relations between them are only weakly dependent on their internal 
structures; any gas equation of state leads to essentially the same relations. We 
briefly consider the possibility that molecular clouds can be described by polytropic 
pressure-density relations (of either positive or negative index), but show that these 
are inconsistent with the apparent gravitational virial equilibrium, 2 U + W ~ 0, 
of GMCs and of massive cores. This class of models would include clouds whose 
nonthermal support comes entirely from Alfven wave pressure. The simplest model 
consistent with all the salient features of GMCs and cores is a "pure logotrope," in 
which P/P c = 1 + A \a{p/ p c ). Detailed comparisons with data are made to estimate 
the value of A, and an excellent fit to the observed dependence of velocity dispersion 
on radius in cores is obtained with A ~ 0.2. 

To appear in the Astrophysical Journal 

1. Introduction 

Giant molecular clouds (GMCs; M ~ 10 5 — 10 6 Mq) in the Galaxy are highly inhomogeneous: 
they are often filamentary in appearance, consisting of discrete clumps, or cores (M <J 1O 3 M0), 
which contain most of the mass of a cloud (including any young stars) and are surrounded by a 
more diffuse component of predominantly atomic gas (e.g., Williams, Blitz, & Stark 1995). Indeed, 
GMCs are clumpy on all scales observed, and are possibly even fractal in nature (Falgarone, 
Phillips, & Walker 1991). Remarkably enough, however, the gross properties of cloud complexes 
are rather simply interrelated. Total masses, mean densities, and average velocity dispersions vary 
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with sizes (effective radii) roughly as M oc i? 2 , p avc oc and cr avc oc R 1 ! 2 (Larson 1981; Sanders, 
Scoville, & Solomon 1985; Solomon et al. 1987), with uncertainties in the exponents typically of 
order ±0.1. The seeming universality of these results demands a physical explanation. 

The relationship between size and linewidth (which term we use interchangeably with velocity 
dispersion) is further interesting because a is observed to decrease towards smaller radii inside 
GMCs, and within individual dense cores (Larson 1981; Miesch & Bally 1994; Fuller & Myers 
1992; Caselli & Myers 1995). Since the total linewidths of GMCs are mostly nonthermal (the 
clouds are stable against gravitational collapse on the largest scales, but have masses several 
orders of magnitude above the thermal Jeans value, so that their support must come largely from 
nonmagnetic or, very likely, MHD turbulence), this decrease reflects a move towards domination 
by thermal motions on the smallest scales. We should therefore expect the linewidths of low-mass 
cores to have a smaller turbulent component than those of high-mass cores. This is indeed the 
case; in fact, the velocities nearest the centers of small cores are almost (to within a few percent) 
wholly thermal. However, it also happens (Caselli & Myers 1995) that the nonthermal velocity 
dispersion shows a stronger dependence on radius in low-mass cores (where cjnt oc r - 5 ) than in 
massive ones (cjnt r 0,2 ). Because star formation is localized in the cores of GMCs, their overall 
structure — and this aspect specifically — bears strongly on our understanding of this most 
fundamental process. 

The goal of this paper is to find a model for the internal structure of molecular cores (low- and 
high-mass both) which quantitatively matches their observed, internal velocity-dispersion profiles, 
and is also consistent with global properties (such as the observed mass-radius-linewidth relations) 
of large, self-gravitating clumps and even whole GMCs. Our approach is to reduce this problem to 
the specification of a (total) pressure-density relation — an equation of state — that, when used 
to solve the equation of hydrostatic equilibrium, results in a gas cloud with the required features. 
It is significant that the linewidth profiles of cores are insensitive to the presence or absence of 
young stars, and thus may be viewed as one of the prerequisites for star formation (Fuller & Myers 
1992; Caselli & Myers 1995). This justifies our focus on the structure of purely gaseous clouds. 

To proceed, we shall resort to a virial-theorem treatment of molecular clouds which idealizes 
them as spheres of gas in hydrostatic equilibrium and satisfying Poisson's equation. This is 
appropriate enough for dense cores, which in many cases are roughly spheroidal (probably prolate: 
Myers et al. 1991), or even near spherical (e.g., Williams et al. 1995), overall. (Although their 
internal density distributions may not be especially smooth on very small scales, we concern 
ourselves here with a description of their bulk structure.) In addition, observations of very massive 
cores imply that they are self-gravitating and in approximate "gravitational virial equilibrium," 
2 U + W ~ 0. This is not the case for low-mass cores, but these still appear to satisfy the full 
virial theorem if surface-pressure terms are included (Bertoldi & McKee 1992). Finally, even 
the total masses, radii, and linewidths of entire GMC complexes are generally consistent with 
virial-equilibrium models of spheres (e.g., Solomon et al. 1987; Elmegreen 1989); and it has been 
repeatedly confirmed that whole clouds tend to comply with 2W + W~0 (Larson 1981; Solomon 
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et al. 1987; Myers & Goodman 1988b). 

It should be noted that the correct, average mass-radius and size-linewidth relations can be 
recovered in a purely scale-free description of GMCs (e.g., Henriksen 1991). But the distinctly 
different scalings of velocity dispersion with radius inside low- and high-mass clumps is one 
indication that molecular clouds are in fact not entirely featureless. Other current models for 
these objects (such as isothermal spheres or negative-index polytropes) face similar difficulties in 
accounting simultaneously for their global properties and their internal structures. 

Magnetic fields are an important presence in regions of star formation (e.g., Heiles et al. 1993), 
so we begin in §2 by writing down the virial-theorem (or Bonnor-Ebert) relations between the 
masses and total linewidths — including turbulent velocities — of magnetized spheres, truncated 
at radii such that the internal pressure just balances that of a diffuse surrounding medium. We 
then develop a stability criterion for such clouds which depends only on the assumption that 
their thermal linewidths (i.e., kinetic temperatures) are invariant. When combined with the 
circumstantial evidence for equipartition between the kinetic and mean-field magnetic energies in 
GMCs and massive cores (Myers & Goodman 1988a, b; Bertoldi & McKee 1992), this stability 
criterion leads to mass-radius-linewidth relations between critical-mass objects that agree with 
the observed scalings (both the exponents and the coefficients) among GMCs. Our analysis is 
therefore similar to those of, e.g., Chieze (1987), Fleck (1988), and especially Elmegreen (1989), 
but ours holds for clouds with an arbitrary equation of state. In addition, we find that critically 
stable clouds in magnetic equipartition should all satisfy 2W + W~0. We therefore conclude that 
GMCs and massive cores are approximately at their critical masses and magnetically supercritical, 
with M crit ^2M$ (see also McKee 1989; Bertoldi & McKee 1992). The generality of these results 
allows for the investigation of essentially any pressure-density relation as a potential description of 
the interiors of molecular cores, just so long as critically stable configurations are at all possible. 

In §3, we discuss polytropic equations of state: P oc p 7 , with 7 any positive number. 
Clearly, if 7 > 1, then the total velocity dispersion (a 2 = P/p) stays constant or decreases with 
decreasing density (increasing radius) inside a cloud. If instead 7 is allowed to be less than 1 
(e.g., Maloney 1988), then a 2 increases with radius, as required. However, in our analysis, any 
such "negative-index" polytrope would be unconditionally stable against gravitational collapse; or 
equivalently (but independently of any stability criterion) , it could not self-consistently satisfy the 
relation 2 U + W ~ 0. While this is not a problem for low- mass cores, it is inconsistent with what 
we know of very large clumps, and GMCs overall. Thus, we argue that polytropic pressure-density 
relations give an incomplete picture of interstellar clouds. A corollary to this is that weakly 
damped Alfven waves, for which P oc p l l 2 (McKee & Zweibel 1995), cannot be invoked as the sole 
explanation of nonthermal linewidths in GMCs. 

Section 4 describes what is, in our view, a more suitable alternative. There we consider the 
possibility that pressure varies only logarithmically with density: P/P c = \ + A ln(p/p c ). We refer 
to the resulting gas cloud as a "pure" logotrope. This term was introduced by Lizano & Shu (1989; 
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see also McKee 1989), who actually added a logarithmic term to an otherwise isothermal equation 
of state in an attempt to account for turbulent linewidths. Such models have also been studied 
in detail by Gehman et al. (1996). Here, however, we dispense with the explicitly isothermal 
component, for two reasons: (1) Assuming P = pa\ + k ln(p/p re f), with <tt the thermal velocity 
dispersion and p re f some reference density (Lizano & Shu 1989), the observational inference that 
linewidths should be essentially purely thermal at the centers of cores requires that p re f = p c . But 
then the internal a 2 = P/p decreases with increasing radius. (2) If now P/P c = p/ p c + n \n(p/p c ), 
then for large values of k such as those suggested by Gehman et al. (1996), P vanishes for p/p c 
rather near unity. Thus, real clouds would have to be almost uniform-density, which again is not 
observed. 

These difficulties do not extend to the specific equation of state which we examine. Instead, 
an outwards-increasing velocity-dispersion profile obtains for an equilibrium pure logotrope. Such 
a model can moreover account for the observed mass-radius and size-linewidth relations between 
GMCs. We also demonstrate that the linewidth measurements in both low- and high- mass cores, 
from a variety of molecular clouds, are quantitatively reproduced if A ~ 0.2. The argument makes 
explicit use of the fact that small cores are not at their critical masses (while large ones generally 
are), but that they are still in approximate virial equilibrium when the effects of surface pressures 
are considered. 

Although the logarithmic P — p relation we advocate is phenomenological, its overall viability, 
along with the failings of other models, ultimately makes for a useful description of molecular 
clouds and their cores. 



2. Generalized Bonnor-Ebert Relations 

Interstellar clouds can be viewed as essentially "pressure-truncated" bodies of gas. (This 
term is meant to imply the existence, not of some radius where a gas cloud suddenly ends, 
but of one where it "blends in" with an ambient medium.) Thus, the boundary of a GMC 
is set by pressure balance with the surrounding, hot ISM. Note that even if there is no such 
balance initially, it will eventually obtain after an overall expansion or contraction of the cloud 
complex. The extent of a core within a GMC is similarly limited by the pressure of a tenuous 
interclump medium, as has been argued by Bertoldi & McKee (1992). The equilibrium structure of 
pressure-truncated isothermal spheres was first described by Ebert (1955) and Bonnor (1956) (see 
also McCrea 1957). Their expressions connecting the masses, radii, and linewidths of the spheres 
follow from the assumption of hydrostatic equilibrium and Poisson's equation, and therefore are 
written in terms of an internal pressure profile and gravitational potential. In Appendix A, we 
re-derive these Bonnor-Ebert relations, but for clouds satisfying an arbitrary gas equation of state 



(eqs. [A4|-[A7]). We also provide a connection (eqs. [ |A8(| -[ A10 |) with the more transparent and 
observationally convenient virial-theorem formulation given by equations ( |2.8| ) below. As various 
authors have noted (Chieze 1987; Fleck 1988; Maloney 1988; Elmegreen 1989), these relations in 
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either guise provide a framework for an understanding of the standard "Larson's laws" for GMCs 
(M oc R 2 and a oc R 1 / 2 ; Larson 1981). 

The development of Appendix A refers specifically to "nonmagnetic" clouds, meaning only 
that no ordered (mean) magnetic field is considered to be present. We work on the assumption 
that the effects of disordered fields (MHD turbulence) can be separately dealt with, in an 
equivalent hydrostatic problem that makes use of an effective equation of state to describe all of 
the contributions to gas pressure as a function of density. Still, the effects of mean fields must also 
be considered in any applications to real interstellar clouds. 



2.1. Magnetic Equilibria 

In the absence of any analytic models for magnetohydrostatic clouds, we proceed by assuming 
spherical symmetry and turning to the scalar virial theorem: 

2 U{\ - Ps/P^e) +M + W = . (2.1) 

Here the mass-averaged, total one-dimensional velocity dispersion of a cloud with radius R is 
related to its mean pressure and density by 



r2 _ Jo fl 47rrV 2 dr _ P a . 
J^4Trr 2 p dr Pa 



so that the kinetic (or internal), mean- field magnetic, and gravitational energies are given by the 
usual 

U = - Mal e , (2.3) 

M = — I B 2 dV + — I(r ■ B)B • dS - — I B 2 r ■ dS , (2.4) 
8ir J An J 8tt J 

and 

„ f R mdm 3 GM 2 , . 

W = -G = - -a ——. 2.5 

jo r 5 R 

The parameter a is essentially a measure of the non-uniformity of a gas sphere, and as such 
depends on the equation of state and the truncation radius (in terms of a fixed scale rg; eq. 
However, it is generally of order unity: for a power-law density profile p oc r~ p , equation (|2.5| ) 
gives a = (1 —p/S)/(l — 2p/5). We expect 1 <^ p <^ 2 in a realistic GMC or core, and thus 
10/9 £ a & 5/3. 

It is often more useful to work in terms of the virial parameter of Bertoldi & McKee (1992): 

a mag = ■ (2.6) 
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This observable quantity can also be written (cf. McKee & Zweibel 1992) as 

a mag - 2a |w| - l _ ^1 |w| 

Although this form of the virial theorem is appropriate for clouds of any shape, the corrective 
factors for spheroidal clouds are rather near unity (Bertoldi & McKee 1992), and we still allow only 
for spherical symmetry here. It is further convenient to distinguish between the virial parameters 
which would obtain for a cloud with and without a mean magnetic field; these are related by 

a mag = CKnon ^1 — jyyj"^ " (2-7) 



Again, a mag refers specifically to the combination ( [2.6D of observables. On the other hand, a non 
applies to the M. = (no mean field) counterpart of a given cloud; it is directly observable only in 
this special case (since then ct mag — c^non); but may always be calculated for a given gas equation 
of state, as outlined in Appendix A. Once this is known, a mag follows with the specification of a 
mean- field configuration (see §2.3). 



The connection ( |2.7D between a mag and a non is valid insofar as the ratio -P s /-P ave does not 
change drastically upon the "addition" of a mean magnetic field to a cloud which is already 
in hydrostatic equilibrium. This must hold along field lines anyway (force balance is required 
in that direction), and thus everywhere on the surface of a roughly spherical cloud. Defining 
a cS = o(l — W|) , such reasoning implies that (a mag — a e ff)/«m ag = (a n on — o)/o:non> so 
manipulation of the virial theorem gives the following: 



3 / Qnnn Q>\ ^ 1 0~ avc 



3 / Qnnn 1 



2 



v /20 ( a non \ 1/2 1 (Ps\ 1/2 



and 



37r \a n(m — a) n l l 2 \G 

«m ag 



Pave 



These relations, of which only two are independent (they are contained, for example, in eqs. [9] 
and [10] of Elmegreen 1989), are also given by Harris & Pudritz (1994) for the specific case of 
critically stable isothermal spheres (as defined in §2.2; a non = 2.054 and a = 1.221 at a cloud 
radius Rcnt/fo = 2.150). As written here, they apply to any generic gas cloud, stable or unstable, 
isothermal or not. 



Any cloud in hydrostatic equilibrium and satisfying Poisson's equation has a central region 
where the potential, density, and velocity dispersion are very nearly constant with radius. 
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Equations flAg ) and (|A9|) show that spheres truncated in or just outside this central region (i.e., 
at radii small enough that P s m P c ) have a non 3> a, and thus a mag esc M~ 2 / 3 oc R~ 2 by equations 
( |2.8| a, b). Roughly this scaling has been observed for the virial parameters of low-mass, high-a 
cores in several GMCs, leading Bertoldi & McKee (1992) to argue that such clumps can be viewed 
as truncated spheroids which are essentially in pressure equilibrium with an intracloud medium. 



2.2. Stability: Critical Clouds 

Given equations fl2.8| ) for the equilibrium structure of pressure-truncated gas spheres, we are 
in a position to question their stability: Under what conditions will they be able to withstand the 
combined effects of self-gravity and surface pressure, and when will they be unstable to wholesale 
gravitational collapse? The answer to this depends, of course, on any boundary conditions 
attached to a perturbation of the cloud. Obviously, the total mass should be unchanged by a 
contraction or expansion of the entire structure. In addition, following Maloney (1988), we suppose 
that the central velocity dispersion remains constant as the cloud radius, or the surface pressure, is 
varied. This stipulation is meant to reflect the fact that the turbulent linewidth decreases steadily 
towards smaller scales in cores and in entire GMCs. We therefore identify the central velocity 
dispersion with the thermal part of the total linewidth: cr 2 = kT/pLniH- Insisting that this be 
invariant amounts to recognizing the rough uniformity of kinetic temperatures T ~ 10K (which to 
first order can be understood as a consequence of the competition between cosmic-ray heating and 
CO cooling) over a large range of scales in interstellar clouds. The stability criterion that follows 
ultimately leads to a set of results which self-consistently explain some important observational 
features of GMCs and massive cores. 

A gas cloud will be stable against radial perturbations if the derivative dP s /dR (taken with 
a c and M held fixed) is < 0: a slight decrease in the cloud radius then leads to an increase in the 
pressure just inside its boundary, which in turn leads to reexpansion. Appendix B shows that, for 
any equation of state, this condition is just 



o- 



9R J m,o> R 



1 - (5/6)(a non - a) 1 

3 Pave/ Ps 



< , (2.9) 



where p s is the internal density at the edge of the cloud. Although this stability criterion has been 
derived without explicitly considering the effects of mean magnetic fields, we expect that it should 
not be greatly altered by their inclusion. (This is again implied by our assumption of approximate 



spherical symmetry, since eq. [ 2^9 1 must at least be satisfied along field lines at the boundary of a 
magnetized cloud.) 

Depending on the equation of state, there may exist a radius for which a pressure-truncated 
cloud is marginally stable (the expression [p.9| is just 0), and beyond which it is unstable. It is 
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these critical equilibria, which must have 

5 

a n on - a = - , (2.10) 

D 

that are of particular interest here. Once a non and a are known as functions of radius (as in 
Appendix A), the satisfaction, if possible, of equation ( |2.10| ) sets the boundary R cv n/rQ of the 
cloud. This in turn allows for evaluation of the coefficients in equations (|2.8|). Quite generally, 



tnont-imag 

1 al 



Rcvit ~ 1/2 1/2 (np\l/2 ' 

V»7TQ!nonQ;mag K^^s) 1 



J crit 
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and 



6ce non P s , ■ . 
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C7. 



avc 



Equation ( p.ll| c) shows that a nonmagnetic cloud (a mag = a non ) on the verge of gravitational 
collapse has a mean column density which is fixed by the pressure of the surrounding medium, 
independently of any gas equation of state. 

In general, the virial parameter of a given cloud is sensitive to its internal structure (through 
the equation of state) and its total radius. However, for a critically stable cloud we always have 
«non = a + 5/6. Since a is typically of order (but slightly greater than) unity, this implies a non ~ 2 
and 1 — Pg/Pgyg = a/a non ~ 1/2. The virial theorem (|2.1| ) then becomes 



U + M + W ^0 . (2.12) 

If there is equipartition A4 ~ IA between the magnetic and kinematic energies in such a cloud, 
then 2 U + W pa as well; thus, A4/|W| ~ 1/2, and the critical a mag is expected to be of order 
unity (cf. eq. [2.7|; see also Elmegreen 1989 and McKee & Zweibel 1992). 



Observations of molecular clouds show that 2 U + W ~ (or equivalently, a mag 1), 
and are consistent with ~ IA (Myers & Goodman 1988a, b), although actual magnetic- field 
measurements are few and uncertain. The most massive cores in GMCs similarly tend to show 
2 U + W ~ and a mag near 1 (this is not the case for low-mass cores, however: e.g., Williams, 
de Geus, & Blitz 1994; Williams et al. 1995; see also §4.1 below). Observations of them are also 
indicative of magnetic equipartition (Bertoldi & McKee 1992), though again the evidence is rather 
indirect, and not necessarily conclusive. Having said this, it does seem that GMCs and massive 



cores both satisfy equation (2.12), which would imply that they are at or near their critical masses. 



Thus arrived at, this conclusion depends on the criterion one adopts for cloud stability, i.e., 
it follows from the assumption that the thermal linewidth a c is held fixed during any radial 
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perturbation (and from the additional proviso that A4 ~ IA). Nevertheless, our analysis — which 
specifies no equation of state — provides a natural explanation for the fact that so many GMCs 
and large cores appear to be in simple "gravitational virial equilibrium" (i.e., 2W + W~0), even 
though they are threaded by appreciably strong mean magnetic fields. 

There is also separate evidence for the criticality of massive molecular cores. For instance, 
Bertoldi & McKee (1992), in their study of the cores in four GMCs, argue that the most massive 
are at least magnetically supercritical, a necessary condition for gravitational instability. More 
fundamentally, star-forming regions clearly must be susceptible to gravitational collapse; but in 
the Rosette GMC at least, those cores which are most obviously associated with IRAS sources are 
also among the heaviest (Williams et al. 1995). As mentioned above, the largest cores in several 
GMCs also have the smallest virial parameters (as low as 1), so that this would seem to be a 
feature of clouds which are close to instability. The observation of a mean a mag « 1 for GMCs 
then implies that they, too, are near some critical mass, and certainly in excess of the nonmagnetic 
Jeans or Bonnor-Ebert value (for which M. = implies IA + W ~ and cc ma g = cn non ~ 2). McKee 
(1989) has further argued that GMCs on the whole are magnetically supercritical, and of course 
they must be strongly self-gravitating in order to be molecular at all (e.g., Elmegreen 1985). 

On a related note, McKee (1989) points out that GMCs must generally be near criticality 
because they show a P ave which is typically an order of magnitude larger than the total (thermal 
plus turbulent) pressure in the hot ISM. Given our stability criterion, equation ( [2,ll| d) shows that 
self-gravity can supply a maximal pressure enhancement P ave /P s ~ 2.5 of a spherical cloud over 
its surrounding medium, and this only for a critical-mass body. Even putting the nonsphericity of 
GMCs aside, however, it is important to note that the analysis here speaks only to the molecular, 
self-gravitating parts of GMCs, and not to their diffuse, low- Ay H I components. Elmegreen 
(1989) has shown that the weight of these atomic "envelopes" can easily increase the pressure at 
the boundaries of the molecular parts of a cloud complex by a factor of 5 or more above the value 
in the ISM at large; overall, then, -P av e/-PiSM > 10. 

If GMCs and their most massive cores are indeed critical-mass objects, then they must all 
have the same dimensionless radii (although the physical scale ro will generally vary), and the 
same virial parameters. Aside from possible variations in P s , which are discussed in detail by 
Elmegreen (1989), this causes the coefficients in equations ( |2.11 ) or (|2.8| ) to be roughly constant, 



and allows for well defined mass-radius-linewidth relations between clouds. Moreover, we have 
argued that a non — 2 and, in the event of magnetic equipartition, a mag ~ 1 for critical clouds, 
regardless of the underlying equation of state. The coefficients in the M — R — <7 avc scalings are 
then independent of this detail, and the average properties of GMCs can shed no light on their 
internal structure. This both explains why critical, magnetized isothermal-sphere models are 
successful in quantitatively accounting for the observed scalings (with a non = 2.054 and a mag ~ 1: 
Elmegreen 1989; Harris & Pudritz 1994), and implies that the same agreement with the data 
comes with any model which provides for the existence of a critical mass (see eqs. 4.12| below). 
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Finally, equations ( |2.11| ) as written would suggest that all GMCs must be under similar 
surface pressures P s if, for example, S is to be roughly the same among them. Other authors 
(e.g., McKee 1989; Mouschovias 1987; Myers & Goodman 1988b) have argued instead that either 
£ itself (and hence -Pave); or the mean field strength B ave , is the more fundamentally invariant 
attribute of GMCs. If so, then P s could be eliminated from the critical Bonnor-Ebert relations in 
favor of any of these quantities, and our approach does not preclude the others. 



2.3. Magnetic Field Model 

Consider a cloud of radius R, threaded by a mean magnetic field approximated as uniform 
and of magnitude -B ave . Outside the cloud, let B fall off as r~ 2 to a radius Rq, where it matches 
onto an ambient, uniform field of strength Bq (Nakano 1984). Conservation of flux (or continuity 
of the normal component of B across the boundary of the cloud) demands B ave R 2 = BqRq, and 
evaluating equation (|2.4| ) at the surface r = R$ gives 

M = B ™ cR * (l- — 
3 V Ro 

Defining j3 = 87rP ave / B 2 ve and $ = ttB^R 2 , we have 

M 5 $ 2 / _ R\ 2a^ f _ R\ 

\W\ 9avr 2 GM 2 \ Ro) 3a/J ^ Rq) ' V ' 

so that equation fl2.7D gives 

1 1 2 / R\ , 

+ • ( 2 - 14 ) 



«mag «non 3a/3 \ Rq J 

With B ~ 3 ^G and S ave ~ 30 - 40 ^G for GMCs (Myers & Goodman 1988b), flux conservation 
gives R/Rq c± 0.3; and U ~ M. implies ~ 1. Thus, a mag rj 1 when a n0 n ~ 2, as expected. 



Some indication of the reliability of equations ( [2.11| ) and ( 2.13 ), and of the approximations 



leading to them, can be had by comparing the critical masses they predict for magnetized, 
isothermal spheres (a = 1.221) with those obtained from self-consistent, axisymmetric numerical 
calculations. In particular, the mass M$ which separates magnetically sub- and supercritical 
clouds, and for which A4 = |W|, is given by equation pl3|) as 0.18 Q/G 1 / 2 , only a 50% 



overestimate of the exact result M$ ~ 0.12 (fr/G 1 / 2 (Mouschovias &: Spitzer 1976; Tomisaka, 
Ikeuchi, & Nakamura 1988). Further, for M £ 0.24 $/G 1/2 , the critical masses we obtain by 
using equation ( p.7| ) in ( |2.11| a) lie within a factor 2 of those found by Tomisaka et al. (1988; 
see their eq. [4.7]). (In fact, our formula is more accurate in the weak-field limit because it 
approaches the correct M cr i t = 1.182 cj 4 /(G 3 P s ) 1 / 2 for the nonmagnetic, = isothermal sphere.) 
And for the equipartition (3 = 1 seen in GMCs, setting a mag ~ 1 in equation ( p. 13 ) implies 



Merit ~ (h/&n 2 ) l l 2 <& /G 1 / 2 . This is roughly 1.6 times our (approximate) M$, and 2.4 times the 
exact value, which level of agreement is quite acceptable. In any case, we are led to expect that 
critical-mass GMCs and cores are strongly magnetically supercritical, with M cr i t ~ 2 M$ — a 
result which has also been argued by McKee (1989) and Bertoldi & McKee (1992). 
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3. Polytropic Equations of State 

The outwards increase of linewidth both within giant molecular clouds as a whole (Larson 
1981; Miesch & Bally 1994) and within individual dense cores (Fuller & Myers 1992; Caselli 
& Myers 1995), immediately suggests the class of negative-index polytropes as possible models 
for these structures. That is, if P oc p 1+1 / N , then for a polytropic index N < —1 we have P/p 
increasing for decreasing p, as required. Such models have been studied by, e.g., Viala & Horedt 
(1974) and Maloney (1988), and we derive the generalized Bonnor-Ebert relations for them in 
Appendix C. Although there are various physical arguments to support their use (e.g., Shu et 
al. 1972; de Jong, Dalgarno, & Boland 1980; McKee & Zweibel 1995), these polytropes turn out 
to be unlikely descriptions of real GMCs (and of course, a positive polytropic index is undesirable 
because it is inconsistent with a velocity dispersion that increases with radius). 

Here we define n = N/(N + 1), so that P oc p 1 /™ and n > 1 for N < —1. From Appendix C, a 
nonmagnetic, negative-index polytrope which is truncated anywhere outside of its constant-density 
central region (inside which, o ~ 1 and a non increases without bound towards r = 0) will then 
satisfy 

5 

a n on - a = -(An - 3) . 

O 

Thus, an n > 1 polytrope has (5/6)(a non — a) -1 < 1; and since its density profile is p oc r~ p , with 
p < 2 everywhere for any n (Appendix C), we also find p me / ' p s < 3. According to equation fl2.9|) , 
then, a truncated, negative-index polytrope will always be stable, and never critically so (as was 
also noted by Maloney 1988). Ultimately, the same steady increase of linewidth with radius which 
would recommend the polytropic equation of state in the first place also proves to be its undoing: 
the internal pressure gradient which results is so shallow as to stabilize a cloud under any external 
pressure. The concept of a critical mass is then irrelevant for these models, which casts doubt on 
their utility in describing real GMCs or high-mass cores. Any additional support from a mean 
magnetic field in the cloud obviously serves only to exacerbate this problem. 

Again, this result follows to some extent from the constraint that the thermal linewidth a c be 
fixed during a perturbation of the cloud. By contrast, both Viala & Horedt (1974) and Chieze 
(1987) consider the possibility that the constant of proportionality in the relation P oc p 1 /™ is 
invariant, and find that polytropes can become unstable for certain truncation radii. Still, our 
approach is closely related to a point which is independent of any rule for cloud stability: 

The minimum value of the nonmagnetic virial parameter for a truncated polytrope is (eq. |C7| ) 

5 (4n-3)(2ra- 1) 
anon " 2 6^5 ' 
As usual, this is a lower limit because a non can be very large indeed if the polytrope is truncated 
at a very small radius. In the event of equipartition U = Ai between kinematic and magnetic 
energies, equation ( |2.7| ) and the identity a mag = 2a U/\W\ lead to 

(V — — 

mag 1 + a non /(2a) 1 + (1/2)(1 - Ps/P^Y 1 ' 
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in which case equation ( pll ) and n > 1 finally imply that 

(4n-3)(2n-l) 10 
amag " (6n-5)(6n + l) 7 ' 

Thus, regardless of whether there exist any unstable modes for truncated polytropes, the virial 
parameters of such clouds are significantly larger than the a mag — 1 (2 U + W ~ 0) seen in GMCs 
(Myers & Goodman 1988b) and in the most massive molecular cores (Bertoldi & McKee 1992; 
Williams et al. 1994, 1995), even if dynamically significant mean magnetic fields are allowed to 
be present. It is a property of our specific stability criterion that this fact implies the absence of 
critically stable equilibria. 

The virial parameters of magnetized polytropes could be reduced to a mag — 1, for any n, if 
the mean field were such that Ai ~ 2 U. In this case, however, using equation QC11 ) for P s /P ave 



in the virial theorem (2.1) leads to 0.8 < .M/|W| < 1 for n > 2, which is difficult to reconcile with 
the rather higher degree of magnetic supercriticality that is observationally inferred for GMCs and 
massive cores (M/\W\ £ 0.5, and M ~ 2M$; see §§2.2, 2.3). 

One consequence of all of this is that the nonthermal linewidths in GMCs, 
°nt = °ave ~~ kT/pmn, cannot be attributed entirely to the pressure of weakly damped 
Alfven waves, for which P oc p 1 / 2 (n = 2, or JV = —2; McKee & Zweibel 1995), and thus 
a mag i; 1-65 under magnetic equipartition. It seems almost certain that Alfven waves do play a 
significant role in the support of GMCs and cores (e.g., Arons & Max 1975; Pudritz 1990); but, 
as McKee & Zweibel (1995) also note, they cannot be uniquely responsible for their large-scale 
stability. In this context, we note that the size-linewidth relation between clouds can be expressed 
in terms of a mean magnetic field strength, as in, e.g., Myers & Goodman (1988a). Specifically, 
equation ( |2.11| d) can be used to write ( |2.11| b) in terms of P ave rather than P s , and the definition 
of f3 (§2.3) relates P ave to -B ave - Then, with a mean mass per particle p = 2.33 and a kinetic 
temperature T = 10 K, we have 

; , i , i I _ ^ave i ( \ 

l magA 



cj nt ~ 0.60 km s" 1 (a mag /3) 1/4 



,30 fiGJ Upc, 

The equality is not quite exact here, because the scalings with l? avc and R strictly apply to 
the total <7 ave . Still, this relation is accurate in the typical case, o"^ T S> kT/pmn] and for 
«mag = /? = 1, it is consistent with available data (see Myers & Goodman 1988a). Thus, although 
Mouschovias & Psaltis (1995) argue for an interpretation of this result in terms of Alfven waves, 
we see here that it is independent of any assumptions on the physical origin of the nonthermal 
motions in GMCs (in principle, they need not even derive from magnetic fields). 

Finally, it is clear that any simple "mixing" of two polytropes with different indices will still 
preclude the existence of a critical (or low-a mag ) cloud; even the superposition of an isothermal 
part, i.e., P = C\p + C2p 1//n , only admits one if it is essentially isothermal anyway {C\ 3> C2). 
Thus, we now turn to a different equation of state, in which the gas pressure varies only 
logarithmically with density. 
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4. The Logotrope 

Many theoretical models of star-forming clouds employ the singular isothermal sphere, in 
which p oc r~ 2 . However, there is some evidence that H II regions are concentrated towards the 
CO centroids of their parent GMCs such that their three-dimensional number density is most 
consistent with an r _1 fall-off (Waller et al. 1987; Scoville et al. 1987). This suggests the possibility 
that, in a heavily smoothed, average sense, the internal density structure of molecular clouds is 
essentially p oc r _1 (see also Solomon et al. 1987). Further, extinction measurements indicate 
p oc r _1 or so in the outer parts of cores as well (Cernicharo, Bachiller, & Duvert 1985; Stiiwe 
1990). 

Lizano &; Shu (1989; see also Gehman et al. 1996) introduced the so-called logotropic equation 
of state for GMCs, P = Pi so + Pturb> with Pi so oc p and -Pturb oc ln(p/p re f). If the central linewidth 
of a cloud is to be entirely thermal in origin, the reference density in Pturb must be p re f = p c - 
However, P/p then decreases with radius, which is incompatible with the observations. We suggest 
that a more complete description of GMCs and cores is given instead by a "pure" logotrope: 



P = PcVc 



1 + A In I — 



(4.1) 



with A > a parameter to be adjusted. The assumption here is that any nonthermal motions, 
which presumably arise from MHD turbulence, add to the thermal pressure such that both are 
fully accounted for by the equation of state ( j4.1| ). This relation is shown schematically in Fig. [l], 
along with an isothermal-sphere combination and the polytropic P oc p 1 / 2 for Alfven waves. In the 
rest of our discussion, equation ( flip is referred to simply as a logotrope, with the understanding 
that no further consideration is given to any mixed-isothermal version. 



Given equation (4J.), the equation of hydrostatic equilibrium (Al) integrates to 

P _ 1 

Pc 1 + ' 



(4.2) 



where if) = {(f) — cf) c )/a 2 is a dimensionless gravitational potential. Substitution of this expression 
into Poisson's equation (|A2|) shows the existence of a singular solution, 



P_ 

Pc 




(4.3) 



which describes the density profile of the bounded (finite p c ) solution at large radii. (The scale 
radius r$ = 3cx c /[47rG'/9 c ] 1 / 2 .) In spite of such a slowly decreasing internal density, a logotrope has 
only a finite extent, since for any positive A the pressure will eventually vanish. Assuming that A 
is small enough for the singular density profile to apply at the edge of the cloud, we have 

Uax = ^ = M e" A . (4.4) 
r V 9 
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The velocity dispersion, a 2 = P/p, increases with radius inside the cloud until it reaches a 
maximum of 

a 2 ma Ja 2 = Ae 1 ^- 1 (4.5) 

at £ = r/ro = e _1 £ max , beyond which it falls again (to 0, at £ max )- 

A logotrope truncated by the pressure P s of an external medium has precisely one critical 
mode for any given value of A. As discussed in Appendix B, if a cloud is truncated at a very small 
radius £, it will always be stable. But in the power-law part of the logotrope, equations (4.2) and 



flOD show that dijj/d£ = (9 A/2) 1 / 2 , so once £ is so large that P s /P c < A/4, the cloud is unstable 



(cf. eq. |B6|). Thus 



-Rcrit 2A f 1 1\ -1/4 1 , | r 



where we have again assumed that A is relatively small and the singular solution (4.3) holds at 
£ cr it- (This expression shows that a achieves its maximum inside £ C rit> and the surface of the 
critical logotrope is cooler than its interior. McKee [1989] has suggested that such behavior might 
arise from the radiation of Alfven waves into the ambient medium.) Now, evaluating the pressure 
( |4.1| ) and density ( |4.3| ) at the radius ( |4.6| ) yields, for any critically stable logotrope, 

Evidently, self-consistency in the use of equation ( |4.3| ) here requires A <C 4. Equations (A7) and 



(2.2) can then be used to find 



3 / al vc \ 14 /P av( ,\ 7 



2 I a 2 ) .. 9 IP, 3 

\ ° / crit 



(4.8 



independently of A. 



A numerical solution for our logotrope with A = 0.18 (which we justify shortly) yields the 
internal density, pressure, and velocity-dispersion profiles shown in Fig. |2[ Equations ( |4.4[ ) and 
( |4.6| ) are valid here, as an r _1 density profile is realized well before the cloud ends (the vertical 
line in all four panels of the Figure marks the radius £ c rit)- It is the steep pressure gradient near 
Cmax which sets the logotrope apart from the negative-index polytropes and allows for a critical 
cloud with a small virial parameter. 

The total velocity dispersion a 2 in Fig. || is roughly constant near the center of the cloud; and 
its rise with radius further on is consistent with a 2 oc r 2 / 3 , which is essentially the scaling originally 
found for GMCs by Larson (1981). Identifying the central dispersion a c with the thermal part 
of the linewidth (so that thermal motions dominate on the smallest scales, as is observed), the 
nonthermal contribution to cloud support is o"^ T = a 2 — a 2 , which is also shown in Fig. |2[ At 
small to moderate radii, <tnt naturally rises more steeply than the total a, while at larger £ the 
two are more comparable in magnitude. Roughly, then, our model has <r^ T oc r over some range 
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in radius before it flattens to cr^ T oc r 1 / 2 and eventually turns over (the maximum 0n T /<t 2 is just 
a Taax/ a c ~ 1) an d OCCUrS at £ = e~ 3 / 4 £ cr ;t). It is found that <tnt ^ r ^ ' within GMCs (e.g., Myers 
& Goodman 1988b), and low-mass molecular cores show the same general trend; high-mass cores 
are more suggestive of ctnt oc r 0,25 (Fuller & Myers 1992; Caselli & Myers 1995). A logotropic 
equation of state both for individual cores and for entire cloud complexes is consistent with these 
observations, provided that low-mass cores can be viewed as spheres which are pressure-truncated 
at radii £ < £ C rit- We discuss this point further in §4.1 below. As always, though, no matter how a 
varies with r internally, the average relation cr ave oc R l / 2 Pl^ still holds for critical-mass clouds, 
i.e., for GMCs. 

A handle on an appropriate value of the parameter A for GMCs as a whole, can be gotten by 
reinventing them as smoothed-out spheres with a number of distinct overdense regions (the cores) 
sprinkled throughout. The observed contrasts in mean density and velocity dispersion between 
GMCs and cores can then be related to the A of a model logotrope. Harris & Pudritz (1994; see 
their Table 2) list each of p ave > oWe; & n d S for a GMC of mass 

M G mc = 3.3 x 10 5 M and a 
core with M core = 5.4 x 10 2 M@. (Such values are "typical" in the sense that, by mass, half of all 
GMCs [cores] are larger than Mqmc pWcore]-) Defining p, a, and E as the core-to-GMC ratios of 
mean volume densities, linewidths, and column densities, these data show that 

p ~ 240 a ~ 0.29 E ~ 4.64 . (4.9) 

Cores with M = M COTe are among the most massive observed, and are likely near critical. If they 
are also logotropes, then their p aV e/p s , etc., should be roughly given by equation ( [4.8|) , regardless 
of whether ^4gmc = A COTe . Such large cores might further be expected to lie in very dense, cold, 
highly pressured regions of their parent clouds. Then the pressure P s at the surface of a core 
is comparable to the GMCs P c , and similarly for the density and velocity dispersion. These 



1/2 

assumptions, together with the relation E oc P s (eq. [2.11c]), result in 



V Ps / GMC K&s/ GMC \ "s J GMC 



Comparison of equations (|4.9|) and (4.7) then shows that the observed p, a, and E (only two of 



which are actually independent) are realized in a logotropic GMC with 

A GMC - 0.175 . (4.11) 



Although this procedure is highly idealized, it is perhaps telling that an A exists at all 
which gives at once the correct p and a. The logotrope appears to be the simplest barotropic 
equation of state which can account for these ratios and the small virial parameters (or the critical 
equilibrium) of GMCs and massive cores. For example, if GMCs and large cores were both n = 2 
polytropes truncated at £ = r/rQ ~ 21.4, then equation (|4.9[) could be roughly satisfied, but they 
would also have a large a non = 5.37, and hence a mag ~ 1.75 for magnetic equipartition (again, 
A4 — 2 IA is required to reduce a mag to unity for a negative-index polytrope, but this is not wholly 
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satisfactory for other reasons; see §3). Alternatively, adding an explicit p 1 / 2 term to equation fl4.1| ) 
destroys the simultaneous correspondence between p and a for critical clouds. 

To this point, we have not accounted for any mean magnetic fields in GMCs and cores. 
However, these should have no effect on our estimate of ^4gmc insofar as very massive clumps 
have the same a mag ~ 1 as their parent clouds. More important is the assumption that the cores 
here lie at the rare, highest-pressure peaks of cloud complexes; if instead they are found in less 
extreme regions, then equation ( |4. 11 ) is an upper limit to Aqmc- Even so, more direct estimates 



of ^4gmc ar e in fairly good agreement with the representative value quoted above. For example, 
equations (4.7) and (4.8) can be combined to obtain an expression for cr 2 ve /a 2 ; comparison with 



observed values of the total cr ave and the thermal a c in GMCs then suggests 0.12 <J ^4gmc ^ 0.16. 

On a final note, since p oc r _1 at the edges of such logotropes, we have 
a = (1 — l/3)/(l — 2/5) = 10/9. With the ratio of mean kinematic and magnetic pressures (3 ~ 1, 
equations ( 2. 10| ) and ( p.!4| ) then give for the critical virial parameters, 



a non = 35/18 and a mag ~ 1.07. 

These are near 2 and 1, as expected from general arguments, and not far from the values 2.054 
and 1.15 for a Bonnor-Ebert isothermal sphere. As alluded to earlier (§2.2), the mass-radius and 
size-linewidth relations among real GMCs can therefore be understood in terms of this model: 
equations (|2.11 b, c) become 



and 



0.37 km s 1 pc 



R 1 / 2 ' v \PismJ V10 4 fecm- 3 K, 

which agree well with observations (cf. Elmegreen 1989). Here P s is the surface pressure on 
the molecular part of the GMC, while Pjsm is the total pressure of the hot, intercloud medium. 
The ratio of these two is typically around 5 — 10, due to the weight placed on a GMC by its 
UV-shielding atomic layer (Elmegreen 1989; also §2.2 above). 



4.1. Comparison With Data: Molecular Cloud Cores 

As a class, dense molecular cores seem somewhat more heterogeneous than their parent 
GMCs. First, many studies (Carr 1987; Loren 1989; Stutzki & Glisten 1990; Lada, Bally, & Stark 
1991) show that the size-linewidth relation between cores can be significantly weaker than the 
virial, <r ave oc R 1 ! 2 scaling among GMCs, and in some cases is difficult to discern at all. Second, 
in the L1630 GMC at least, the size-linewidth relation is considerably better defined for clumps 
which stand out more strongly against the interclump medium (see the comparison in Fig. 13 of 
Lada et al. between their "5cr" and "3cr" cores). And third, the virial parameters a ma g of GMC 
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cores are not restricted to values near 1, but rather can range as high as ~100 (Bertoldi & McKee 
1992; Williams et al. 1994, 1995). 

This last point is particularly significant because the virial parameters of cores in a number 
of different clouds are seen to correlate with their size: low mass goes along with rather large 
a mag , and high mass with more moderate a mag ~ 1. Thus, as Bertoldi & McKee (1992) argue, the 
smallest clumps in GMCs are far removed from gravitational instability (surface pressure is more 
important than self-gravity in their confinement), while the largest are actually near criticality 
(strongly self-gravitating; see also §2.2 above). Equivalently, a range in the masses of cores can 
be viewed, at least roughly, as a range in the dimensionless radii £ — small for low mass, and 
near £ cr it for high mass — at which they are truncated by the ambient pressure in a GMC (virial 
parameters generally decrease with £ for any equation of state; Appendix A). This notion could 
account for the rather confused size-linewidth relation between cores (if £ does not have a fixed 
value, then neither do a non and a mag in eq. |2.8| b1), and ultimately must also have consequences 
for the interpretation of data on the radial variation of velocity dispersion inside cores. 

If indeed very low-mass cores differ from very high-mass ones mainly in being truncated at 
Clow ^ £high — Ccrit; then the former must be less centrally condensed, have ratios P s /P c nearer 
unity, and be generally less distinguishable from the intercore (GMC) gas than the latter. This 
results in the scale radius ro = 3cr c /(47rG / o c ) 1 / 2 oc ct 2 /P c 1/ ' 2 being larger for less massive clumps. To 
see this, compare the ro of a low-mass core, truncated at £ small enough that P s P c , with the 
ro of a near-critical, high-mass core. On average, any two cores should be under similar surface 
pressures P s (this essentially being set by the internal P ave of a parent GMC), so that 



ro(high) fP s \ 1/2 



r (low) \P C J high 



CT 2 (high) 


[a 


CT 2 (high) 


_ of (low) _ 


V 4 


_ o-c(kw) _ 



(4.13) 



where the second step follows from equation fl4.7|) if a logotropic equation of state applies. 
Thus, for a c not too widely different between the two cores, and for A small enough, the ratio 



(4.13) will be < 1. The same physical radius r in high- and low-mass cores then corresponds 
to respectively larger and smaller £ = t/tq. This result, which just reflects the fact that more 
strongly self-gravitating structures are more centrally concentrated, must be considered when 
attempting to match any model to any observations. 

Returning now to the issue of internal velocity-dispersion profiles, Fuller & Myers (1992; 
hereafter FM) and Casein & Myers (1995; hereafter CM) have compiled linewidth measurements 
in at least three different molecular lines for each of 14 low-mass and 24 massive cores. The 
resulting a vs. r profiles indicate that, over similar ranges r ~ 0.1 — 1 pc in all the clumps, and 
independently of whether or not stars are present within them, the velocity dispersion rises more 
steeply with radius in low-mass cores than in high-mass ones. In particular, nonthermal linewidths 
in the former are consistent with cjnt oc r 0,53 ; in the latter, ctnt oc r ' 21 (CM). In our view, this 
situation is a direct result of the effect summarized by equation ( 4.13| ). That is, we consider the 



low-mass cores of FM to be eminently stable logotropes truncated at small £, while the high-mass 
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clumps of CM are near critical stability and hence have smaller r^. The two surveys target similar 
r within each type of core, so the massive-core data apply to larger dimensionless £, and must 
sample the shallower parts of velocity-dispersion profiles like those in Fig. 

To actually confront the predictions of our logotropic equation of state with the observations 
of FM and CM, we again identify the (model) central velocity dispersion with the (observed) 
thermal linewidth: a 2 = o\ = kT/fimn, where T is taken from CM, and \i = 2.33. CM assume 
T = 10 K for all the low-mass cores, and we further assume them all to have P s ~ P c ; they 
should therefore have a common tq, which may be easily estimated. Five of the FM cores 
are observed at radii ttnt such that cjnt = ct; among these, (txnt) — 0.12 pc (see Table 3 
of CM). The related quantity £tnt = ttntAo in a model logotrope is just that point where 
a 2 /a 2 = (o"^ T + a^)/aj, = 2. If, for example, A = 0.2 in equation ( |4.1[) , we find £tnt = 0.51, and 
thus 

ro(low) ?s 0.25 pc . 

This result is then used in equation ( 4.13| ), along with a 2 oc T and A = 0.2, to give 

' rri \ 

J- high 



ro(high) « 0.056 pc 



10 K 



where 12 K < T < 33 K for the massive cores (CM). The data may now be appropriately scaled 
and compared to theoretical a — r and <tnt — r curves. 

Figure ^ is the result of such a comparison with three logotropes of different A. The good 
overall agreement between the models and data suggests that GMC cores are well described by 

A core = 0.20 ± 0.02 , 

which is encouragingly close to Aqmc as given in equation ( 4.11| ) above. In Table |l| we list the 



radii and other properties (as calculated from eqs. 4^] to [}4.7f above) of critical logotropes with 
A near 0.2. 

Some of the scatter in Fig. || could arise if A does not have a common value for all cores, 
and any embedded stars could further complicate the situation. (Notably, FM and CM both find 
that starless cores alone show a tighter — but not different — linewidth-radius relation than do 
cores with stars, suggesting that the basic form seen here is one of the initial conditions for star 
formation.) However, some of the scatter is surely due to the fact that the cores in this sample 
do not reside in a single GMC, but come from diverse environments. (For instance, the low-mass 
cores here tend to be found in dark clouds of mass ~ 10 4 - 1O 5 M , and the high -mass ones m 
somehat larger GMCs.) Thus, there is no guarantee that every low-mass core observes P s ~ P c , 
or that every high- mass core is at its critical mass. It is the large number of cores available here 
which allows the mean trend in Fig. |||, and the implied structural dichotomy between small and 
large clumps, to emerge. 

As with our earlier consideration of core-to-GMC density and linewidth ratios, a fit of similar 
quality to the data in Fig. can be obtained with the model velocity-dispersion profile of an n = 2 
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polytrope which has a scale ro about a factor of two smaller than that used here. As we have 
stressed, however, there are other difficulties with the application of polytropic equations of state 
to GMCs and cores. 

In summary, the approach we have taken self-consistently allows for purely thermal motions to 
dominate on the smallest scales in cores (and in GMCs); for an internal density structure (p oc r _1 
in the outer parts of cores, and much shallower near their centers) which is consistent with at 
least some observations^ (Cernicharo et al. 1985; Stiiwe 1990; Andre, Ward-Thompson, & Barsony 
1993); and for a unified treatment of low- and high-mass cores. This is in contrast to the models 
of FM and CM (also Myers & Fuller 1992), which assume two distinct components of isothermal 
and nonthermal gas in cores and predict a p oc r~ 2 singularity at their centers. It should also be 
noted that the assumption P s ~ P c in low-mass cores still allows for these noncritical clumps to 
show a fair degree of central concentration. For example, an A = 0.2 logotrope which is truncated 
at £ ~ 0.1 £ cr it has p s /p c < 0.09, but P s /P c ~ 0.5; thanks to the relatively weak dependence of P 
on p, we do not require, nor even expect, that all low-mass cores be uniform-density spheres. Our 
model is therefore consistent with the significantly non-uniform density profiles of low-mass cores 
(a mag £ 3) in the Rosette GMC (Williams et al. 1995). 



5. Summary 

We have generalized the classic Bonnor-Ebert relations for isothermal spheres, to give 
expressions for the masses, radii, and total velocity dispersions of magnetized, pressure-truncated 
gas clouds with any internal pressure-density relation. The analysis combines an exact approach, 
based on solving the equation of hydrostatic equilibrium, with a virial-theorem treatment to 
incorporate mean magnetic fields into the clouds. A stability criterion has been developed that 
relies only on the assumption of an invariant central velocity dispersion a c , and is independent 
of the gas equation of state. Clouds which are just critically stable under this condition have 
mass-radius and size-linewidth relations which are effectively oblivious to their internal structure, 
i.e., to their equation of state. 

These results have been applied to molecular cloud complexes and dense cores, leading to 
three main conclusions: 



1 Williams et al. (1995) fit the projected density profiles of many clumps in the Rosette GMC, with a function of 
the form (1 + r p /a)~" — a being the projected half-power radius. They find n ~ 1 outside of unresolved central 
regions, and suggest that this implies an intrinsic density profile of p oc r~ 2 . However, because molecular cores have 
a finite extent, their density profiles in projection are not so simply related to their space densities. (For example, 
the surface density of a core with any p profile must decline steeply in the outermost regions, as less and less of the 
core material is intercepted by the line of sight.) Even a truncated logotrope, for which we might naively expect a 
flat surface-density profile, actually gives rise to something consistent with n ~ —0.8 to —1.0 over most radii in the 
fitting function of Williams et al.; and the projected half-power radius is only a = 2.92ro ~ 0.2 pc for a critical, 
A = 0.2 core (R crit = 24.37r ). 
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(1) In spite of their nonspherical geometry, the gross features of GMCs are consistent with 
those of critical-mass clouds in our analysis. (In the context of GMCs and cores, we identify the 
central velocity dispersion with the thermal part of the total [thermal plus turbulent] linewidth.) 
Independently of the exact equation of state, if equipartition between magnetic and kinetic energies 
obtains in a critically stable cloud, then it will appear to be in simple gravitational equilibrium, 
2 U + W ~ (equivalently, it will have a virial parameter a mag ~ 1). These are observed features 
of real GMCs and their most massive cores. In addition, any of our critical-mass spheres obey the 
same mass-radius and size-linewidth relations which hold among real GMCs. All of this implies 
that GMCs, and large clumps, are near criticality. 

(2) No polytropic (power-law) pressure-density relation can fully characterize GMCs or cores. 
A positive polytropic index leads to a velocity dispersion which decreases outwards within a cloud, 
opposite to what is seen both in individual cores and in entire cloud complexes. Alternatively, 

a negative index (i.e., P oc p 7 with 7 < 1) gives a linewidth which increases with radius, as 
required. However, such clouds do not show virial parameters near unity unless their mean fields 
are so strong that their magnetic and gravitational energies are comparable in magnitude, and 
this is again in conflict with observations of GMCs and large cores. Given our stability criterion, 
an equivalent statement is that pressure-truncated polytropes are unconditionally stable against 
gravitational collapse. Since weakly damped Alfven waves satisfy P oc p 1 / 2 (McKee & Zweibel 
1995), this shows explicitly that they cannot fully explain the observed nonthermal support in 
molecular clouds. 

(3) The most successful model for the internal structure of molecular cores, and one which is 
also consistent with global properties of GMCs, is a "pure logotrope:" P/P c = 1 + A ln(p/p c ). 
This relation is meant to account for all contributions to the total gas pressure, including the 
effects of disordered magnetic fields (MHD turbulence). Once it is recognized that low-mass 
cores are far below their critical masses (though still in virial equilibrium; they are essentially 
pressure-confined), and that high- mass cores are near criticality, the internal velocity-dispersion 
profiles of clumps from a variety of environments are seen to be consistent with a logotropic 
model with A ~ 0.2. A similar value of A can also explain characteristic core-to-GMC ratios of 
mean densities and linewidths, and any A is consistent with the observed size-linewidth relation 
between GMCs because the logotropic equation of state allows for critical equilibria. Finally, the 
equilibrium density profile of a logotrope has p oc r _1 , outside of a constant-density central region. 
There is some observational support to be found for this prediction, in GMCs and cores both. 

We would like to thank Phil Myers for his comments on an earlier version of this paper, and 
Charles Curry for valuable discussions. R.E.P. also acknowledges conversations with Paola Caselli 
and Gary Fuller. This work was supported in part by the Natural Sciences and Engineering 
Research Council of Canada. 
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A. PRESSURE-TRUNUCATED EQUILIBRIA 

For a spherical cloud in hydrostatic equilibrium, with kinematic (thermal plus turbulent) 
pressure P = pa 2 (so that a is the one-dimensional velocity dispersion) and a self-gravitational 
potential 4>, we define 

i> = - f c , r o = A 9 ^ > and £ = r / r a, 
of. 4:irGp c 

where a subscript c denotes evaluation of a quantity at the cloud center. The characteristic scale 
of a cloud is set by ro; the factor of 9 in its definition identifies it with the projected half-power 
radius of an isothermal sphere (e.g., Binney & Tremaine 1987). The equation of hydrostatic 
equilibrium then reads 

d fP\ p dtp 

and Poisson's equation becomes 



d€\Pj Pc di ' (A1) 



1 d ( 2 dip\ n p 



^wJ-V (A2) 

(Note that ip = dip/dtl = at £ = 0, and tp > for £ > 0.) The mass enclosed within radius £ is 
given by 

M ^ M(£) = 4 Wo 3 /V-<£' , 

j0 Pc 

so that, with the help of equation QA2j) , 

M -T / ' cro l e ^eJ - l e rfeJ • (A3) 

Following Ebert (1955) and Bonnor (1956), the surface of the cloud is defined by that radius 
£ at which the internal pressure P(£) just equals a confining pressure P s due to a surrounding 
medium of negligible gravity. Equation (|A3|), the definition of ro, and the identity P c = p c & 2 then 
yield, for any equation of state relating P and p, 



"=W*Mi?J wfa- (A4) 



The physical radius of the cloud is R = £ro: 



These relations combine to give 



M PL dnj> (P s y l/2 (Ps\ 1/2 ( \r\ 

S = ^ = V97 deUJ UJ (A6) 
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and 

3M 1/1 c#\ fPsY 1 P s 



Given an equation of state, the integration of equations (Al) and (A2) determines dip/d£ and 
P s /P c at any truncation radius £. This fixes the coefficients in equations (|A4| )~([A7|) and affords 
mass-radius and size-linewidth relations for pressure-truncated clouds. (Note that the scalings in 
these relations — M oc a /Ps /2 , R oc a 2 /Ps^ 2 , etc. — are global ones, and do not necessarily 
reflect the radial dependence of any quantity inside a cloud.) 



The definition of the non-uniformity parameter a (eq. |2.5| ) allows it to be found as a function 
of radius: 

15 fdip\' 2 ft A/3 dV P 



'f (A8) 

io a£ p c 



e we. 

and equation (|A3|) implies that the virial parameter 5a 2 ve R/GM for nonmagnetic clouds (i.e., 

C^mag — C^non) IS just 

«-= 5 f («!)"• (A9) 

which tends to be a nonincreasing function of £. Alternatively, equations (|A7| ) and ( |2.8| d) give 



a aon - a = 15 \ d 4r Y ( § ] • (A10) 



These results can be used to rewrite equations (|A4|) -( A7) in terms of a n o n! a, and a ave ; equations 
(|2.8| ) (with a mag = a n on) are then obtained. Also, once a non is known as a function of radius in 
any cloud, the observable a ma g for a given magnetic field configuration may be calculated (e.g., 
§2.3). 



B. STABILITY CRITERION 



We begin by recalling the definition of the scale radius tq from equation (A3): 



GM / 2 dip 
<r£ V d£ 



(Bl) 



and rearranging equation (A4): 



47T 



dip 



Ps = ^U Z1 T7 "FT 



P, 



cr ; 



Pj G 3 M 2 



(B2) 



The stability of a cloud truncated at radius R by the external pressure P s is determined by the 
sign of the derivative (dP s /dR): this is negative for a stable equilibrium, for the "critical" cloud, 
and positive for instability. 



-23- 



We now apply a perturbation R — ► R + d-R, while keeping the mass M and central (i.e., 
thermal) velocity dispersion a c of the cloud fixed, so that 



dPs 
dR 



fdP s 



M,a c 



M,a c 



dR 



M,cr c 



(B3) 



Then R = £ro and equation (p3l|) give us 



or) 



M,ct c 



1 

TO 
1 

ro 
1 

?•() 



^0 V 9^ / A/, CTc 



d£/ d£ 



1-9$ 



\dU Pc 



-I -1 



(B4) 



where p s is the density just inside the cloud edge and we have used Poisson's equation (A2). We 
also have, from equation ( |B2[) , 



(dP s 



Ps 



M,cr c 



df 



d 2 ^ 4 cft/A fP s \ fdi)\ 2 d 



de + i da \p c 



d£J dip VP 



Ps\-rr\ - 



Pc 



dip 



18- ^ 



p 



Pr 



(B5) 



where both equations (Al) and ( |A2| ) have been used. 
Thus, equations dB3|), (p4]), and dB5[) show that 



Up; 



= 18^ 

M,(7 C P 



l-(l/18)(^M) 2 (Ps/P c ) 



-1 



(l/0(dl/>/dt)(p c /ps)-9 



(B6) 



Now, from equations ( |A10| ) and (A7), we have 
1_ (d4\ 2 (PA' 1 = 5 1 



so that finally, 



18 \ d£,J VP 



'dP M 



6 a r 



V 3P 



M,(T C 



- 6 i 



i /dv 



1 - (5/6)(a non - a) 1 



3 Pave/ Ps 



(B7) 



Alternatively, for a nonmagnetic cloud, we can make use of equation ( |2.8| c), with a mag 
write 



^noni to 



( — ) 

\dRj 
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M,cr c 



P 



1 - (7rGS 2 )/(8P, 

3 — Peeve/ Ps 



(B8) 



This treatment is valid for any gas equation of state, so long as the cloud radius is perturbed in 
such a way that M and o~ c are unaffected. The result flB8| ) was also obtained by Ebert (1955) and 
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Bonnor (1956), who specifically considered the stability of truncated isothermal spheres, and by 
Maloney (1988) in his study of negative-index polytropes. 

Any virialized cloud truncated at small £ (i.e., where p and o 2 are approximately constant) 
will be stable. This is because, regardless of the equation of state, the potential ip at small radii 
can always be expanded in a power series of the form: ip = (3/2)£ 2 — 0(£ 4 ) (see Appendix C for an 
example). Thus, referring to equation (B6) with p s p c , we have (dP s /dR)M,a c ~ —3P S /R < 0. 
Depending on the equation of state, there may or may not be a truncation at larger £ which 
results in (dP s /dR)M,a c > an d an unstable equilibrium. 



C. POLYTROPES OF NEGATIVE INDEX 

Consider the (nonmagnetic) equation of state 

P<xp 1+ V N , 

where the polytropic index N < — 1. Then we define n = N/(N + 1) and write 

P = p c a 2 c (P-) 1/n , n>l, (CI) 



\pc/ 

with n = 1 corresponding to isothermality. In order to evaluate the coefficients in equations 
flA|-([A7|) (or £|, with 

O'mag — c^non)) we first use the relation ( pl[ ) to integrate the equation of 



hydrostatic equilibrium (Al) and find 



£. = [i + ( n _ (C2) 

Pc 



so that Poisson's equation (|A2| ) becomes 



l±f e W\ =Q[1 + in _ lW /(l-n) (C3) 



i 2 d£ V d£ , 

(Note that in the limit n — * 1, the right-hand side of eq. |C3(| becomes 9e - ^, just what is required 
for the isothermal sphere.) At the cloud center, if) = dip / dtl = 0, so for small £ the potential ip can 
be expanded in a power series, 



allowing for an integration of equation ( p3[ ). The result of this is a core-halo structure to the 
cloud: at very small radii, the solution is indistinguishable from that for an isothermal sphere, and 
the density and velocity dispersion are essentially constant with £, independently of n. Once £ ^ 
a few, however, a power-law density profile obtains; for each n, there exists a singular solution to 
which the equilibrium cloud tends at large radii. 

These singular solutions are readily found by substituting the prescription 

JL = KC P 

Pc 
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into equation ( p2j ) and solving (C3) for K and p. This gives 



p = and K 

1 2n - 1 



H (P-1)(3-P) 



p/2 

(C5) 



so that with n > 1 we have 1 < p < 2, and these clouds are diffuse enough that they have no 
"natural" edge. That is, the density never vanishes, and any boundary to the cloud must be 
defined by the point at which the internal pressure matches that of an external medium. (It is 
again worth noting the exact correspondence to the singular isothermal sphere when n = 1 in 



eq. [C5|; cf. Chandrasekhar 1967.) 



Consider now a polytrope of negative index, which is truncated at radius £ large enough 
that the structure of the cloud is given by the singular solution. Then dtp/dS, and P s /P c can be 
calculated analytically, and equation ( A1C|) yields 



5 

a n on - a = -{An - 3) . (C6) 



We also know that a = (1 — p/3)/(l — 2p/5), so 



6n 

and we have the following: 



5 (4n-3)(2n-l) 

"non = 7T a ? V^>) 



M -Vn{An-3r (L^l) {G^K) 1 ' 2 ' (C8) 



V 27r (4n-3) 2n-l (GP a )Va ' 



V vr(4n -3) \G 
and 

6n-3 P, 



- 1 , (cio) 



6n — 5 <7 2 . 



(Cll) 



avc 



Thus, regardless of where they are truncated by the surface pressure P s , these model clouds can 
only ever represent modest enhancements over the intercloud medium. In particular, P ave /P s has 
a maximum of 3 (for isothermality, n = 1), but even for n = 2 (which gives the p dependence of 
Alfven wave pressure; McKee & Zweibel 1995) is reduced to just 9/7. 

As discussed in §3, any truncation of a negative-index polytrope results in a stable cloud. 
This is because, with p oc r~ p , we have that p a ve/Ps = 3/(3 — p), and the stability criterion ( |B7|) 
becomes (with the help of |Cq| ) {dP s /dK)M,a c = —AP S /R < 0. In particular, this holds in the limit 
n — * 1, when the numerator and denominator of ( [B7D both vanish. What allows for the existence 
of critical and unstable equilibria for the bounded isothermal sphere, then, is the fact that its 
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density profile actually oscillates about the singular solution p oc r~ 2 at large radii (Chandrasekhar 
1967). This does not occur for those polytropes with n > 1; instead, the bounded spheres follow 
the singular solutions exactly at large £. 

Finally, it is worth noting that these results do not contradict the well known instability (e.g., 
Shu 1977) of truly singular spheres, for which p oc r~ p all the way to the center and p c — ► oo. In 
such cases, perturbation of the truncated cloud radius is performed with M, a c , and p c all held 
fixed. Stability then depends only on the sign of (dP s /d£)jtf i(Tc (eq. |jB5| ), which is always positive 
for these polytropes. 
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Table 1. Properties of Critical Logotropes. 



A 


£max 


£crit 


Ps/pc 






Ps/Pc 


0.18 


51.73 


40.29 


4.96 x 10~ 3 


9.07 


8.07 


0.045 


0.20 


31.29 


24.37 


8.65 x 10~ 3 


5.78 


4.78 


0.050 


0.22 


20.83 


16.22 


1.36 x 10" 2 


4.04 


3.04 


0.055 
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Fig. 1. — Comparison of the logotropic equation of state [P/P c = 1 + 0.18 ln(p/p c )] with those 
for Alfven wave pressure [P/P c = (p/Pc) 1 ^ 2 ] and a logotrope-plus-isothermal sphere combination 
[P/P c = p/p c + 0.1 ln(p/p c )]. 

Fig. 2. — Structure of a logotropic gas sphere with A = 0.18. The vertical line in all four panels is 
at the truncation radius which makes for a critically stable cloud: £ cr i t = 40.29. The long-dashed 
line in the plot of p vs. r is the singular density profile of eq. (4.3); that in the plot of total linewidth 
(bottom left) represents a 2 oc r 2//3 ; and those in the plot of nonthermal linewidth (bottom right) 
trace a^ T oc r and a 2 NT oc r 1 / 2 . Recall that r 2 = 9cr 2 /4irGp c . 

Fig. 3. — Comparison of observed internal velocity-dispersion profiles (total a, thermal and 
nonthermal components o"t and cjnt) with logotropic models, for 38 GMC cores. Open squares 
correspond to low- mass cores (Fuller Sz Myers 1992), and filled symbols to high- mass cores (Caselli & 
Myers 1995; we omit four HCO + measurements). Cores both with and without stars are represented 
here, and each has been observed in three or more different molecular lines. The model curves are 
for A = 0.20 (best case; solid line), and A = 0.18, A = 0.22 (dashed lines); a larger A results 
in smaller maximum and critical radii r/ro- The vertical line is at £ cr it = 24.37 for an A = 0.2 
logotrope. 
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